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Abstract. The perturbation theory is developed based on small parameters 
which naturally appear in solid state quantum computation. We report the 
^ , simulations of the dynamics of quantum logic operations with a large number 

OO 1 of qubits (up to 1000). A nuclear spin chain is considered in which selec- 

CO ' tive excitations of spins are provided by having a uniform gradient of the 

C 3 ' external magnetic field. Quantum logic operations are utilized by applying 

^^ , resonant electromagnetic pulses. The spins interact with their nearest neigh- 

C _ 3 , bors. We simulate the creation of the long-distance entanglement between 

remote qubits in the spin chain. Our method enables us to minimize un- 
wanted non-resonant effects in a controlled way. The method we use cannot 
simulate complicated quantum logic (a quantum computer is required to do 
this), but it can be useful to test the experimental performance of simple 
quantum logic operations. We show that: (a) the probability distribution of 
unwanted states has a "band" structure, (b) the directions of spins in typical 
S ^ unwanted states are highly correlated, and (c) many of the unwanted states 

are high-energy states of a quantum computer (a spin chain). Our approach 
can be applied to simple quantum logic gates and fragments of quantum 
algorithms involving a large number of qubits. 
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j^ ■ Introduction 

Recently much progress has been made in single particle technologies. These 
technologies allow one to manipulate a single electron, a single atom, and a single 
ion. A distinguishing feature of all these technologies is that they are "quantum" ; 
quantum effects are crucial for the preparation of the initial state, for performing 
useful operations, and for reading out the final state. One of the future quantum 
technologies is quantum computation. In a quantum computer, the information 
is loaded in a register of quantum bits - "qubits" . A qubit is a quantum object 
(generalized spin 1/2) which can occupy two quantum states, |0) and |1), and an 
arbitrary superposition of these states: ^ = Co |0) -I- Ci 1 1) . (The only constraint is: 
|CoP + |Cip — 1.) A quantum computer is remarkably efficient in executing newly 
invented quantum algorithms (the quantum Fourier transform, Shor's algorithm 
for prime factorization, Grover's algorithm for data base searching, and others). 
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Using these algorithms, quantum computing promises to solve problems that are 
intractable on digital computers. The main advantage of quantum computation is 
the rapid, parallel execution of quantum logic operations. One of the promising 
directions in quantum computation is solid state quantum computation. (See the 
review pi.) In solid state computers, a qubit can be represented by a single nuclear 
spin 1/2 |2|, 1^, or a single electron spin 1/2 [Q, |5|, a single Cooper pair ||^, 0]. 
This type of quantum computers is quite different from quantum computers in 
which a qubit is represented by an ensemble of spins 1/2 ||^, ^ p!o| . 

Crucial mathematical problems must be solved in order to understand the dy- 
namical aspects of quantum computation. One of these problems is the creation 
of the dynamical theory of quantum computation - the main subject of our paper. 
The processes of the creation of quantum data bases, the storage and searching 
of quantum information, the implementation of quantum logic gates, and all the 
steps involved in quantum computation are dynamical processes. When qubits 
representing a register in a quantum computer are in superpositional states, they 
are not eigenstates of the Hamiltonian describing the quantum computer. These 
superpositions are time-dependent. Understanding their dynamics is very im- 
portant. To design a working quantum computer, simulations of quantum logic 
operations and fragments of quantum computation are essential. The results of 
these simulations will enable engineers to optimize and test quantum computers. 

There are two main obstacles to perform useful quantum logic operations with 
large number of qubits on a digital (classical) computer: (1) the related Hilbert 
space is extremely large, Dn = 2^, where N is the number of qubits (spins), and 
(2) even if the initial state of a quantum computer does not involve many basic 
states (eigenstates) , the number of excited eigenstates can rapidly grow during the 
process of performing quantum logic operations. Generally both of these obstacles 
exist. At the same time, many useful quantum operations can be simulated on a 
digital computer even if the number of qubits is quite large, say 1000. How one 
can do this? One way is to create (and use) a perturbation theory of quantum 
computation. To do this, one should consider a quantum computer as a many- 
particle quantum system and introduce small parameters. Usually, in all physical 
problems there exist small parameters which allow one to simplify the problem to 
find approximate solutions. The objective in this approach is to build a solution 
in a controlled way. 

As an example, it is useful to consider an electron in a hydrogen atom inter- 
acting with a laser field. The electron has an infinite (formally) number of discrete 
levels plus a continuous part of its energy spectrum. So, a single electron in a 
Hydrogen atom has a Hilbert space larger than all finite qubit quantum comput- 
ers. Assume that the electron is populated initially on some energy level(s). The 
action of the laser field on the electron leads in many cases to a regular and con- 
trolled dynamics of the electron. In particular, if the amplitude of the laser field 
is small enough, the electron interacting with this laser field can be considered as 
a two-level quantum system interacting with an external time-periodic field. Why 
is the electron excited to only few energy levels? This happens because the ex- 
istence of small parameters makes the probability of unwanted events (excitation 
of the electron to most energy levels including the continuum) negligibly small. 
The existence of small parameters allows one to use perturbation approaches when 
calculating the quantum dynamics of the electron. This strategy can be used to 
create a perturbation theory of quantum computation. 



In this paper we present our results, based on the existence of small parame- 
ters, for simulating quantum logic operations with a large number (up to 1000) of 
qubits. Our perturbation approach does not provide a substitution for computing 
quantum algorithms on a quantum computer. One will need quantum computer to 
do the required calculations. Our method will allow one to simulate the required 
benchmarks, quantum logic test operations, and fragments of quantum algorithms. 
Our method will enable engineers to optimize a working quantum computer. 

In section 1, we formulate our model and the Hamiltonian of a quantum 
computer based on a one-dimensional nuclear spin chain. Spins interact through 
nearest-neighbor Ising interactions. All spins are placed in a magnetic field with 
the uniform gradient. This field gradient enables the selective excitation of spins. 
Quantum logic operations are provided by applying the required resonant electro- 
magnetic pulses. In section 2, we present the dynamical equations of motion in the 
interaction representation. As well, we discuss an alternative description based on 
the transformation to the rotating frame. In the latter case the effective Hamil- 
tonian is independent of time and solution of the problem can be obtained using 
the eigenstates of this Hamiltonian. In section 3, we derive simplified dynamical 
equations taking into account only resonant and near-resonant transitions (and 
neglecting non-resonant transitions). We introduce the small parameters which 
characterize the probability of generation of unwanted states in the result of near- 
resonant transitions. We describe a 27rA:-method which allows us to minimize the 
influence of unwanted near-resonant effects which can destroy quantum logic op- 
erations. In section 4, we describe a quantum Control-Not gate for remote qubits 
in a quantum computer with a large number of qubits. Analytical solution for 
this gate is presented in section 5, for the 2Trk condition. In section 6, we present 
small parameters of the problem, in explicit form. The equation for the total 
probability of generation of unwanted states, including the states generated in the 
result of near-resonant and non-resonant transition, is derived in section 7. In 
section 8, we present results of numerical simulations of the quantum Control-Not 
gate for remote qubits in quantum computers with 200 and 1000 qubits when the 
probability of the near-resonant transitions is relatively large (when the conditions 
of 27rfc-method are not satisfied). We show that: (a) the probability distribution 
of unwanted states has a "band" structure, (b) the directions of spins in typical 
unwanted states are highly correlated, and (c) many of the unwanted states are 
high-energy states of the quantum computer (a spin chain). The total probability 
of error (including the states generated in the result of near-resonant and non- 
resonant transitions) is computed numerically in section 9 for small {N = 10) and 
large {N = 1000) number of qubits. A range of parameters is found in which the 
probability of error does not exceed a definite threshold. We test our perturba- 
tive approach using exact numerical solution when the number of qubits is small 
(N = 10). In section 10, we show how the problem of quantum computation (even 
with large number of qubits) can be formulated classically in terms of interacting 
one-dimensional oscillators. In the conclusion, we summarize our results. 

1. Formulation of the model 

The mathematical model of a quantum computer used in this paper is based on 
a one-dimensional Ising nuclear-spin system - a chain of N identical nuclear spins 
(qubits). Application of Ising spin systems for quantum computations was first 



suggested in Ref.[ll]. Today, these systems are used in liquid nuclear magnetic 



resonance (NMR) quantum computation with small number of qubits 1 12 . The 
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Figure 1. Nuclear spin quantum computer (the ground state of 
nuclear spins). B^{z) is the slightly non-uniform magnetic field; 
bj^ is the amplitude of radio-frequency field of the n-th rf pulse. 
The chain of spins makes an angle 9 with the direction of the 
magnetic field B^{z). 



register (a ID chain of N identical nuclear spins) is placed in a magnetic field, 
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where t'^^^ < t < i("+i) and n = f, ...,M. In (1), B^{z) is a slightly non-uniform 
magnetic field oriented in the positive z-dircction. The quantities bj^ > 0, v^'^' > 
and if^"^' are, respectively, the amplitude, the frequency, and the initial phase of the 
circular polarized (in the x — y plane) magnetic field. This magnetic field has the 
form of rectangular pulses of the length (time duration) t^"^ = t("+i) — t("). The 
total number of pulses which is required to perform a given quantum computation 
(protocol) is M . Schematically, our quantum computer is shown in Fig. 1. 

1.1. The quantum computer Hamiltonian. The quantum Hamiltonian 
which describes our quantum computer has the form. 



(2) 



E 

k=0 



N-2 



H = Ho + V^-J2^kr,-2jJ2 Ikik+i + V. 



k=0 



(We set h = 1.) The operator V describes the interaction of spins with pulses of 
the r/ field, 



(3) y = ^l/(" 



M 

) 



n=l 



and y*-"' describes the interaction of spins with n-th pulse of the rf field, 






fe=0 



+ 



/,+ 



exp 



(v^'^h + ifi^'A 



where 8*^"-'(i) equals 1 only during the nth pulse. The operators in (||)-(^ have 
the following explicit form in the z-representation (the representation in which the 
operators i| are diagonal), 

(5) il^\\^k){Ok\-\\lk){lkl 

it = H + ^il = lOfcXUI, 4" - it - ^il - luxo^l- 

We shall use the Dirac notation for the complete set of eigenstates (the stationary 
states) of the quantum computer described by the Hamiltonian Hq in (0). The 
eigenstates of the spin chain can be described as a combination of 2^ individual 
states of nuclear spins, 

(6) I00...00), |00...01),...,|ll...ll), 

where the state |0) corresponds to the direction of a nuclear spin along the direction 
of the magnetic field, B^{z), and the state |1) corresponds to the opposite direction. 
The subscript "fc" is omitted. 

In (y)-(y), ujk is the Larmor frequency of the fc-th spin (neglecting interactions 
between spins), LUk = ^B^{zk), and 7 is the nuclear gyromagnetic ratio. (For 
example, for a proton in the field B^{zk) = lOT, one has the nuclear magnetic 
resonance (NMR) frequency fk — ujk/^n w 430MHz.) We assume, for definiteness, 
that the gradient of the magnetic field is positive, dB^{z)/dz > 0. Suppose that 
the frequency difference of two neighboring spins is, Sf = 6uj/2tt « lOkHz, where 
Suj ~ ujk+i — ujk- Thus, if the frequency of the edge spin is 430MHz, the frequency 
of the spin at the other end of the chain of 1000 qubits is « 440MHz. Then, the 
value of B^(z) increases by AB^ = 0.23T along the spin chain. Taking the distance 
between the neighboring spins, a ~ 2A, we obtain the characteristic value for the 
gradient of the magnetic field, dB^{z)/dz « O.23/lOOOacos0, where 9 is the angle 
between the direction of the chain and the z-axis (Fig. 1). Below we will take 
cos 6 = l/\/3- (This allows one to suppress the dipolc interaction between spins in 
the eigenstates of the Hamiltonian Hq.) Thus, the gradient of the magnetic field 
is dB^{z)/dz w 2 X lO^T/m. The quantity fJ^") = "fo"^' is the Rabi frequency of 
the n-th pulse. 
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2. Quantum dynamics of the computer 

We discuss below the dynamics of the spin chain described by the Hamiltonian 
(||)-(||). The quantum state of the quantum computer is described by the wave 
function ^(<). The dynamics of this function is given by the solution of the 
Schrodinger equation, 

(7) i^(t) = il*(t), 

where dot means derivative over time. This linear equation and an initial condi- 
tion, Vl/(0), define the state of the system at the time t. In this section we describe 
two different approaches which allow to solve Eq. ffl). 

2.1. Equations of motion in the interaction representation. In order 
to compute the dynamics it is convenient to go over to interaction representation. 
The wave function, ^i„t(i), in the interaction representation, is connected to the 
wave function, ^(i), in the laboratory system of coordinates by the transformation, 

(8) *(f) = C/o*»nt(t), Uo = exp(-iiJoi), 



where, Hq is defined in (0). We choose the eigenfunctions, \p), and the eigenvalues, 
Ep, of the Hamiltonian Hq as the basis states and expand the wave function ^mt(t) 
in a series, 

(9) ^^nt{t) = E ^pWIp)' 

where the states |p) satisfy the equation 

(10) Ho\p)^Ep\p), p^l,...,2''. 

The wave function in the laboratory system of coordinates has the form, 

(11) *W = E ^pWIp) (=M~^Ept), 

p 

Using the Schrodinger equation (0) for the wave function, $(t), we obtain the 
equations of motion for the amplitudes Cp (t) , 

2"-l 

(12) tCp = J2 ^'""^'p™ e'^P [*(^P - ^™)^ + *^P™ (^'"^^ + ^^"^)] ' 

which represent the wave function in the interaction representation. Here rpm = 
=F1 for Ep > Em and Ep < Em, respectively, Vpm = —^^"^/2 for the states \p) 
and \m) which are connected by a single-spin transition, and Vpm = for all other 

states. 

2.2. The dynamics in the rotating frame. Sometimes, to calculate the 
quantum dynamics generated by the Hamiltonian (0) during the action of the n-th 
pulse, t*-"-* <t< i(""'"^\ it is convenient to make a transformation to the system 
of coordinates which rotates with the frequency j^^"^ of the magnetic field, b± . To 
do this, one can use the unitary transformation |lJ 

(13) [/(") = exp 



7V-1 



^U-h + ^("^)J2U 



k=0 



In the rotating system of coordinates the new wave function, "^rot (t) and the new 
Hamiltonian, Ti'") have the form, 

(14) ^,,t{t) ^ ;/(")+*(*), 



N-l 



(15) 7i:(") = ;/(")+//(") ;/(") = -J2 



k=0 



- j/("M k + f7(")/f + 2jhi' 



LOk-iy^-~' ik+^^' '^k+^-J^k-^k+i 



where the term I^_^I^ must be excluded. In the rotating system of coordinates 
the effective magnetic field which acts on the fc-th spin during the action of the 
n-th pulse has the components, 

(16) i?(") = (f^("),0,c.fc-i.("))/7. 

The advantage of the Hamiltonian (|l^) is that it is time- independent. In this 
case one can use the eigenstates of the Hamiltonian Ti^") to calculate the dynamics 
during the nth pulse without solution of the system of differential equations (112) . 
As before, we expand the wave function in the rotating frame, ^rot{t)^ over the 
eigenstates of the Hamiltonian H^) , 

(17) *.o,(i) = Y. Mt)\p)- 

The wave function in the laboratory system of coordinates is, 

(18) *(i) = Y. ^p(^)Ip) e^p (-*4"^^ + 4"^' 

p 

where x^"' = " (-^"V2) E^"o"' '-f , &^ = W^'V^)i:t'.' ^l < = "1 if the 
state |p) contains the fcth spin in in the state |1), and ct^ = 1 if the state |p) 

contains the fcth spin in the state |0). Below we take (/s*-"' = ^p = for all n. 
The Schrodinger equation for the amplitudes Ap(i\ has the form, 

(19) ^A,[i) = (e, - x(")) ^p(t) - ^- ^ A,,(i), 

v' 

where the sum is taken over the states |p') connected by a single-spin transition 
with the state |p). 

vv' 
effective Hamiltonian, Tt!^"^ is independent of time. When the number of spins, 

vv' 



Equation (|19|) can be written in the form iAp(£) = Tip Mp'(t), where the 

; number of sp 
iV, is not very large (iV < 30), one can diagonalize the Hamiltonian matrix H}^ 



and find its (time-independent) eigenfunctions, 

(20) V^")=^^^("))b) 

V 

and eigenvalues, e>q . The amplitudes ^p "' are the eigenfunctions of the Hamil- 
tonian Ti^") in the representation of the Hamiltonian iJg. 

Using the eigenstates of the Hamiltonian Ti^"-' one can calculate the dynamics 
generated by the nth pulse as, 

(21) ^p(i„) = A,^Y. ^Po E ^?i"^^?^"^ exp (-ze(")r„) , 

Po 9 
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where t„ — tn-i + t„, r„ is the duration of the pulse, Apg — Ap(tn^i) is the 
amphtude before the action of the nth pulse, Ap — Ap(tn) is the amplitude after 
the action of the nth pulse. 

The amplitudes, Ap(t), in the rotating frame are related to the amplitudes, 
Cp{t), in the interaction representation by (see Eqs. ( pi] ) and (p^)), 

(22) Ap{t) = exp(-z4")t)Cp(i), 

Where Ep — Ep — Xp ^-re the diagonal elements of the Hamiltonian matrix Ti-Jpi ■ 
Since the rf pulses in the protocol are different, the effective time-independent 
Hamiltonians, Hll, , are also different for different n. Hence, before each nth pulse 
we should make the transformation to the rotating frame, Cp(t„_i) — + Ap{tn-i), 
and after the pulse we must return to the interaction representation j4p(f„) -—^ 

L'p(in)- 

In the rotating frame the most difficult problem is the diagonalization of the 
sparse symmetric matrices (for each pulse) of the size 2^® 2^. (In each row of this 
matrix there are only 7V-|- 1 nonzero matrix elements.) This approach can be used 
directly (without perturbative consideration) only for small quantum computers, 
with the number of qubits TV < 30. For our purposes, to model quantum dynamics 
in a quantum computer with large number of qubits (TV = 1000), we shall use a 
perturbation approach based on the existence of small parameters. 

3. Resonant and near-resonant transitions 

The number of equations in (|lj) grows exponentially as N increases. Thus, 
this system of equations cannot be solved directly when TV is large. Our intention 



is to solve Eqs. (12) or (|l9| ) approximately, but in a controlled way. This can be 
done if we make use of small parameters which exist for the system. The explicit 
expressions for small parameters will be given below, in sections H and H. 

3.1. Resonant and non-resonant frequencies. In this section, we de- 
scribe the main simplification procedure which is based on the separation of reso- 
nant and non-resonant interactions. This procedure is commonly used in studying 
linear and non-linear dynamical systems |14], 15, 16]. The main idea of this ap- 



proach is the following. When the frequency of the external field, v^^'^ is equal (or 
close) to the eigenfrequency of the system under consideration, Wpm = Ep — E^, 
and when a small parameter exists (the interaction is relatively weak), one can 
separate the slow (resonant) dynamics from the fast (non-resonant) dynamics. The 
main contribution to the evolution of the system is associated with resonant effects. 
Usually, non-resonant effects are small, and can be considered using perturbation 
approaches. In quantum computation, the problem of non-resonant effects is more 
complicated because these effects can accumulate in time creating significant er- 
rors in the process of quantum computation. So, in quantum computation, the 
non-resonant effects must not only be taken into consideration, but they must 
also be minimized. Below, we describe methods which allow us to minimize the 
influence of non-resonant effects. 

In classical dynamical systems, the separation of resonant and non-resonant 
effects is very efficient when "action-angle" variables are used. In this case, the 
"action" is a slow variable and the "angle" (the "phase") is a fast variable. From 
this point of view, Eqs. (O) are convenient because they are written in the "en- 
ergy" representation in which the slow (resonant) effects and fast (non-resonant) 
effects can be naturally separated. 



First, we note that as the number of spins, N^ increases, the number of eigen- 
states, \p), increases exponentially (as 2^), but the number of resonant frequencies 
in our system grows only linearly in N. Indeed, the number of resonant frequencies 
is 3A^ — 2 because only single-spin transitions are allowed by the operator V in 
(^ . The definition of the resonant frequency is the following. We introduce the 
frequency of the single-spin transition: LUp„i = Ep — Em, where Ep and Em are 
two eigenstates, \p) and \m), of the Hamiltonian Hq with opposite orientations of 
only the fc-th spin, (fc — 0, ..., N — 1). For example, 

(23) \p) = \nN^inN-2--0k--ninQ), \m) = \nN-inN-2--lk--nino), 

where Uk = 0,1. The resonant frequency for transition between pth and mth 
states, Vres, of the ?i-th external r/ pulse is defined as follows: 

(24) 4fs = Kra\. 

The only resonant frequencies of our nuclear spin quantum computer are, 

(25) Wfc ± J, (k = 0, or k ^ N - 1), 

Wfc, u;k±2J, {l<k<N~-2), 

where the Larmor frequencies, cuk = jB''{zk)- AH resonant frequencies in ( pq ) are 
assumed to be positive. For end spins with fc = and k = N — 1, the upper and 
lower signs in (|2^) correspond to the states |0) or |1) of their only neighboring 
spins. For inner spins with 1 < A: < A'^ — 2 the frequency ujk corresponds to having 
nearest neighbors whose spins are in opposite directions to each other. The "-f" 
sign corresponds to having the nearest neighbors in their ground state. The "-" 
sign corresponds to having the nearest neighbors in their excited state. 

As an example, consider the resonant frequency of a single-spin transition 
between the following two eigenstates, 

\p) = |nAr_inAr_2---l/c_iO/clA;+i...nino), 

\m) = |nAr_in7v-2---lfc_ilfclfe+i...nino). 

These two eigenstates differ by a transition of the fc-th spin from its ground state, 
\0k), to its excited state, |lfc). All other spins in the quantum computer remain in 
their initial states. According to (0), the eigenvalues of the states \p) and \m) are: 

(26) Ep = E„res -Y+J^ Em^ E^res + ^~J^ 

where E^res is the total energy of all A^ — 1 (except the fc-th) which did not 
participate in the transition. It follows from (|2l) that in this case, the resonance 
frequency, VrVs, of the n-th pulse is, 

(27) v';!ll^\Ep-Em\=0Jk-'2J. 

This frequency belongs to the set of the resonant frequencies introduced in (p5|). In 
the following, we assume that our protocol includes only the resonant frequencies 
from (18), and we will omit the subscript "res". 
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3.2. Dynamical equations for resonant and near-resonant transi- 
tions. In this sub-section we derive the approximate equations for resonant and 
near-resonant transitions, in the system. Consider an arbitrary eigenstate of the 
Hamiltonian Hq, 

(28) |nAr_inAr_2...nino), 

where, as before, the subscript indicates the position of the spin in a chain, and 
rzfe = 0, 1. Assume that one applies to the spin chain a resonant rf pulse of a 
frequency, i^^"\ from (Eq). Then, one has two possibilities: 

1) The frequency of the pulse, i^("\ is the resonant frequency of the state (P8|). 

2) The frequency, i^*^"\ differs from the closest resonant frequency of the state (|2S 
by the value A — ±2 J or ±4J (near-resonant frequency). 

In the first case, one has a resonant transition. In the second case, one has a 
near-resonant transition. If the small parameters exist, J/6lu <C 1, and fl/Sto <C 1, 
we can neglect all non-resonant transitions related to a flip of fc'th non-resonant 
spin. Below, in section O, we will write a rigorous condition for Su which is required 
in order to neglect all non-resonant transitions. 

Thus, considering the dynamics of any state (E8[) under the action of an rf 
pulse with the frequency, i/'^"\ from (p4|), (P5|), we need only to take into account 
one transition. This transition will be either a resonant one or a near-resonant 
one. This allows us to reduce the system of equations (|lj) to the set of only two 
coupled equations (two-level approximation), 

(29) iCp = -(r!(")/2) exp[z(^p - E^ - iy^'^^)t]Cm, 

zC„ = -(17("V2) exp[i(^™ -Ep + i^("))t]Cp, 
where Ep > E„i, \p) and \m) are any two eigenstates which are connected by 
a single-spin transition and whose energies differ by z^*^"^ or z/'"' -I- Apm, where 
Apm = Ep — Em — J^'"' is the frequency detuning for the transition between the 
states \p) and |m). 

3.3. Solutions of the dynamical equations in two-level approxima- 
tion. The solution of Eqs. (P9| ) for the case when the system is initially in the 
eigenstate |m), is, 

(30) 

Cm{to + t) = [cos(Ap„jr/2) + i{Apm/Xprn) sin(Aj„„r/2)] x exp(-irAp„j/2), 

Cp{to + t) = i{Q/Xpjn) sm{XpmT/2) x exp(iioAp„j -I- irAp„/2), 
Cm{ta) ~ 1, ^3(^0) = Oj 

where Ap„i = (fJ^ + A^^)^/^, Iq is the time of the beginning of the pulse; r is 
the duration of the pulse. Here and below we omit the upper index "(n)" which 
indicates the number of the rf pulse. If the system is initially in the upper state, 
b), (C™(to) = 0, Cp(io) = 1), the solution of (||) is, 

(31) Cp(io + t) = [cos(Ap„r/2) - i(Ap„/Ap,„) sin(ApmT/2)] x exp(iTAp™/2), 
Cm (to + t) = i(f7/Ap„) sin(Ap,„r/2) x exp(-itoApm - irAp„/2), 
Cp{U) = 1, C„,(to) = 0. 

For the resonant transition (Ap^, = 0), the expressions ( po|) and (|3^) transform 
into the well-known equations for the Rabi transitions. For example, we have from 
H) for Ap„. = 0, 

Cm(io + r) = cos(rjT/2), Cp(to + r) = isin(f7T/2). 

10 



In particular, for ilr = tt (the so-called 7r-pulse), the above expressions describe 
the complete transition from the state |to) to the state \p). 

For the near-resonant transitions, consider two characteristic parameters in 
expressions (|3^) and (pl|): 

If either of these two parameters, tpm or epJi, is zero, the probability of a non- 

(2) 

resonant transition vanishes. The second parameter, tpm, is equal zero when, 

(32) Xp,nT = 2TTk, (fc=l,2,..), 

where k is the number of revolutions of the non-resonant (average) spin about the 
effective field in the rotating frame. Eq. (p3) is the condition for the "27rfc" -method 



to eliminate the near-resonant transitions. (See |13 , Chapter 22, and |17, UM.) 



3.4. The 27rfc-method. We shall present here the explicit conditions for the 
"27rfc" rotation from the state \p) to the state \m). For a 7r-pulse (fir = tt), the 
values of Q which satisfy the 27rfc condition are, according to Eq. (p2), 

(33) 

nf'i = |Ap„|/V4fc2 - 1 = \Ap^\/V3, |Ap„|/%/l5, |Ap„|/%/35, |Ap™|/V63, ... 
For a 7r/2-pulse, the corresponding values of fi are. 



(34) f^g^l = |Ap„|/Vl6fc2-l = \Apm\/Vl5, \Ap^\/V63, . . . 

(If the Rabi frequency, fl, satisfies the 27rfc condition (|34| ) for a 7r/2-pulse, it 
automatically satisfies the condition ( p3| ) for a 7r-pulse.) 

3.5. Resonant and near-resonant transitions in the rotating frame. 

We consider here the structure of the effective time-independent Hamiltonian ma- 
trix Tippi in the rotating frame. We discuss the two-level approximation in the 
rotating system of coordinates and demonstrate equivalence of the solution in the 
rotating frame with the solution (Bfl) (or (|3l|)) in the interaction representation. 

Let us discuss the structure of the matrix Tipp' . Since the spins in the chain 
are identical, all nonzero non-diagonal matrix elements are the same and equal to 
(— rj/2). At n <^ 6uj the absolute values of the diagonal elements in general case 
are much larger than the absolute values of the off-diagonal elements. We explain 
below how the resonance is coded in the structure of the diagonal elements of the 
Hamiltonian matrix, Tipp'. 

Suppose that the kth spin in the chain has resonant or near-resonant NMR 
frequency. The energy separation between pth and ?nth diagonal elements of the 
matrix Tippi related by the flip of resonant (or near-resonant) fcth spin is much 
less than the energy separation between the pth diagonal elements and diagonal 
elements related to other states, which differ from the state \p) by a flip of a non- 
resonant fc'th {k' y^ k) spin. In this case one can neglect the interaction of the pth 
state with all states except for the state \m)^ and the Hamiltonian matrix Jippi 
breaks up into 2^/2, approximately independent 2x2 matrices, 



(35) 



S-m V 

V £p, 



where Ep = £"„, + Ap„j, |Apm| ~ J or zero, and V = —Q/2 is the perturbation 
amplitude. 



Suppose, for example, that N — 5 and third spin (fc — 3) has resonant or near- 
resonant frequency. (We start enumeration from the right as shown in (p3).) Then, 
the block 2x2 will be organized, for example, by the following states: JOIOIO) 
and 100010); |01111) and |00111); |00001) and |01001), and so on. In order to find 
the state |m), which form 2x2 block with a definite state \p), one should flip the 
resonant spin of the state \p). In other words, positions of iV — 1 (non-resonant) 
spins of these states are equivalent, while position of the resonant spin is different. 

We now obtain the solution in the two-level approximation. The dynamics is 
given by Eq. (^ . Since we deal only with a single 2x2 block of the matrix Tipp^ , 
(but not with the whole matrix), the dynamics in this approximation is generated 
only by the eigenstates of one block. In order to demonstrate equivalence of two 
descriptions (in the rotating frame and in the interaction representation) we will 



choose the same initial condition as in (30). Then, we will make the transformation 



22|) to the rotating frame, Cp —> Ap. After that we shall compute the dynamics by 
Eq. ( pi|) and return to the interaction representation, Ap -^ Cp, using Eq. (|2^). 
Eventually, we will obtain Eq. (po|). The result is the same in the rotating frame 
and in the interaction representation. Physically these two approaches are equiv- 
alent (exactly, but not only in the two- level approximation), but mathematically 
they are different. In the interaction representation we calculate the dynamics 
generated by the time-dependent Hamiltonian. There are no stationary states in 
this case and one should solve the system of differential equations. On the other 
hand, in the rotating frame the effective Hamiltonian Ti is time-independent and 
the wave function evolves in time because it is not the eigenfunction of Ti. 

The eigenvalues Cq , Cq , and the eigenfunctions of the 2x2 matrix ( p5| ) are 
(we put Apm = A and Xpm = A), 

(36) 40)=£„ '^ ' ^^™"'^ ' ( ^ 



2 2' l^4'"7 V(A-A)2 + f}2 \\-A 

(371 e(°)-£ +^ + ^ 1^-4^'°'^ 1 r-(A-A) 

^^^^ '« -f™+ 2 + 2' l^^pQwJ - yiA^AF + TI^ I " 

Suppose that before the nth pulse the system is in the state |m), i.e. the conditions 

C™(io) = 1, Cp(to) - 0, 

are satisfied. After the transformation, (p2), to the rotating frame wc obtain 

Am{U) = exp(-if,„io)C'm(io) = exp(-z£,„to), ^p(io) = 0. 

The dynamics is given by Eq. (|2l|), which in our case takes the form: 



A^(t) = A„,(io) 



A^(°)) exp(-*e(")r) + (^^(")) exp (-ze[,°V 



exp[-i[£-„,t-(A/2)r]] 



{f}2e-.^^/2 + (;^„A)2g.Ar/2|^ 



r!2-f (A-A)2 
where t = to + t- Applying the back transformation, 

C„(i) ^ exp(if„it)A„(t), 

and taking the real and imaginary parts of the expression in curl brackets we 
obtain the first equation (pQ). For another amplitude we have, 

Ap(i) = An(to) [^^(''^A^^") exp (-ze^r) -f A^(")^Q(°) exp (-^egV)" = 



i-exp{~i[£mt- (A/2)T]}sin(AT/2). 
A 

Applying the back transformation. 

Cp{t) = exp[2(f™ + A)t]Ap{t), 

we obtain the second equation in (pw. 

One may demonstrate the equivalence of our two approaches in a different, 
more simple, way. The dynamical equations with the 2x2 Hamiltonian ( pq ) are, 

lAp = tpAp + V Ajn, 
lAm = V Ap + tmAm- 

After the transformation (E2|) to the interaction representation we obtain, 

(38) iCp = Vexp[i{£p - £,n)t]Cjn, 

iCm = Vexp[i{£„i - £p)t]Cp, 

where £p — £m = Ep — Em ~ v — Apm, V — —17/2. One can see that Eqs. ( JSq ) 
are equivalent to Eqs. ( P9|) with the solution given by (JS^) or (pl|). 

4. A quantum Control-Not gate for remote qubits 

The quantum Control-Not (CNki) gate is a unitary operator which transforms 
the eigenstate, 

(39) \nN-i--nk n;....no), 

into the state, 

(40) \nN^i:-nk ni....no), 

where n; = 1 — n; if n^ = 1; and n; — ni if Uk = 0. The fc-th and l-th qubits are 
called the control and the target qubits of the CNki gate. One can also introduce 
a modified quantum CN gate which performs the same transformation (59), (En) 
accompanied by phase shifts which are different for different eigenstates [|13[|. It 
is well-known that the quantum CN gate can produce an entangled state of two 
qubits, which cannot be represented as a product of the wave functions of the 
individual qubits. 

In this section, we shall consider an implementation of the quantum CN gate 
in the Ising spin chain with the left end spin as the control qubit and the right 
end spin as the target qubit, i.e. CA^jv-i.o for a spin chain with large number 
of qubits, N {N = 200 or N = 1000). Using this quantum gate we will create 
entanglement between the end qubits in the spin chain. 

Assume that initially a quantum computer is in its ground state, 

(41) * = |0...0). 

Then we apply a 7r/2-pulse with the resonant frequency, 

(42) i^^LUN-i- 

The 7r/2-pulse means that the duration of this pulse is r = 7r/2ri. This pulse 
produces a superpositional state of the {N — l)-th (left) qubit, 

(43) *^(l/\/2)(|0...0>+i|1...0)). 

To implement a modified quantum CNn^i^q gate we apply to the spin chain 
M — 2N — 3 TT-pulses which transform the state |10 . . . 00) to the state |10 . . . 01) 
by the following scheme: |1000...0) -^ |1100...0) -> |1110...0) -^ |1010...0) -^ 
I1011...0) -^ I1001...0) ->•••-> I100...11) -> I100...01). All frequencies of 

13 



our protocol are resonant for these transitions. If we apply the same protocol 
to the system in the ground state, then with large probability the system will 
remain in the ground state because these pulses have the detunings from resonant 
transitions, A„ = AJ,^ ^ 0. Here the ground state |0) = |00. . .00) is related to 
the state \m) by the flip of the /cth spin with the near- resonant NMR frequency. 
The first 7r-pulse has the frequency lo = a;7v-2- For the second 7r-pulse uj = wat-s. 
For the third 7r-pulse, lo — lun-2 — 2 J, etc. All detunings, A„, in our protocol are 
the same, A„ = 2 J, except for the third pulse, where A3 = 4J. 



5. An analytic solution with application of a 27rfc-niethod 

Assume that all M 7r-pulses satisfy the 27r/c condition, with the same value of 
k: 

(44) A[flT„ = 2Trk, (and n^^^T„ - tt). 



where aI^' = ^(l^lf')^ + A^ In this case, we can derive an analytic expression 

for the wave function, ^, after the action of a 7r/2-pulse and M tt- pulses. This 
solution has the form, 

(45) vE' = Co|00..0) + Ci|10...1), 
where, 

(-l)kM f-l)^-l 

(46) Co = -5^ — ^ exp(--i7rMv/4fc2-l/2), Ci = -^ '-^, 

V 2 V 2 

5.1. Large k asymptotics. For fc 3> 1, we get the same solution for odd 
and even k: 

Co » I/V2. 

This result is easy to understand. For a 7r-pulse, the Rabi frequency is, 

(47) rfW = |A„|/V4fc2-1. 

Large values of k correspond to small values of the parameter: 

^L'V|A„i«i. 

As f2„ /|A„| approaches zero, the non- resonant pulse becomes unable to change 
the quantum state. 

5.2. Small k behavior. For small values of k, the non-resonant pulse can 
change the phase of the initial state. For example, for k — 1 we have, 

(48) Co = exp[i7rM(l - \/3/2)]/\/2. 

After the first 7r-pulse, the phase shift is approximately 24°. This phase shift grows 
as the number of 7r-pulses, M, increases. This increasing phase is an effect which 
can be easily controlled. 
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6. Small parameters 

We shall now introduce the small parameters of the problem. Consider the 
probability of non-resonant transition. This probability will be small if £„ is small: 

(49) £„ = {^n/K? sin2(A„r„/2) « 1. 

(It follows from (|3^) that the expression for |Cmp can be written in the form: 
\Cm\^ = 1 — £„, and \Cp\^ — En-) In order to minimize the errors caused by the 
near-resonant transitions, it is reasonable to keep the values of £„ the same and 
small. Since the values of the detuning are the same for all pulses, A„ = A = 2 J 
(except for the third pulse, where A3 = 4J), and because e„ depends only on 
the ratio |r2„/A„| we take the values of il„ to be the same, r2„ — Q, [n ^ 2i) and 
f23 — 2ri. In this case e„ is independent of n (e„ — e). If we take into consideration 
the change of phase of the generated unwanted state, this change will be small if, 

(50) r!„/|A„|«r!„/A„«l. 

(For the 27rfc-method, this condition requires fc ::^ 1.) 

Next, we will discuss the probability of the near-resonant transitions, after the 
action of M pulses, using small parameters e„. Analytic expressions for probabil- 
ities |Cop and |Cip are: 

M 

(51) |Co|'=2n(l-^"), |Cil'=2- 

n=l 

If all values of e„ are the same for all pulses, £„ — e, then in the first non- vanishing 
approximation of the perturbation theory we obtain, 

(52) \Co\' = lil-Me), |CiP = i. 

The decrease of the probability, |CoP, is caused by the generation of unwanted 
states. One can see that the deviation from the value, |CoP = 1/2, grows as 
the number of 7r-pulses, M , increases. It means that the small parameter of the 
problem is Me rather than e. When M « 2N is large, even a small deviation from 
the 27rfc condition can produce large distortions from the desired wave function 

(il)- 

7. Non-resonant transitions 

In this section we estimate the errors caused by non-resonant transitions and 
write the formula for the total probability of errors including into consideration 
both near-resonant and non-resonant transitions. In order to estimate the prob- 
ability of non-resonant transitions it is convenient to work in the rotating frame, 
where the effective Hamiltonian is independent of time. 

Consider a transition between the states \m) and |m') related by a flip of a 
non-resonant fc'th spin. The absolute value of the difference between the mth and 
TTi'th diagonal elements of the matrix Hpj,' is of order or greater than Scu, because 
they belong to different 2x2 blocks ( p5| ) . Since the absolute values of the matrix 
elements which connect the different blocks are small, \V\ <C 6uj, we can use the 



standard perturbation theory [19|. The wave function, Vjg, in ( pOf ) can be written 
in the form, 

q' S« ~ ^q' 
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where superscript "(n)" is omitted, prime in the sum means that the term with 
q' = q is omitted, V'g is the eigenfunction of the Hamiltonian Ti, the gth eigenstate 
is related to the mth diagonal element and the q'th eigenstate is related to the m'th 
diagonal element, Vqgi — 2V{ipq l^fc/IV'g' ) is the matrix element for transition be- 
tween the states ipq and ip, ' , the sum over q' takes into consideration all possible 
non-resonant one-spin- flip transitions fronr the state Jto). Because the matrix Tipp' 
is divided into 2^~^ relatively independent 2x2 blocks, the energy, eq (eL ), 

and the wave function, Vq {wq' )i in ( |53| ) are, respectively, the eigenvalue ( p^ ) 
and the eigenfunction, 

(54) <) = A^(«)|m) + A^(°)|p), {^pf = ^?:(«)|m') + < ("^b')) 

of the single 2x2 block (35) with all other elements being equal to zero. 

When the block (35) is related to the near-resonant transition (A ~ J) and 
when ^ <^ J ^ the eigenfunctions of this block are, 

(55) ^f « [1 - (f}V32j2)]|m) + (17/4J)|p), 

^g"«-(f^/4J)|m) + [l-(f7V32j2)]b). 

On the other hand, if this block is related to the resonant transition (A = 0), we 
have, 

(56) <) = (l/N/2)(H-fb)), V^g') = (l/V2)(|m)-|p)). 

The probability of the non-resonant transition from the state \vi) to the state |m') 
connected by a flip of the non-resonant fc'th spin is, 

(57) Pm^' = \{'m'\%)f . 

Note, that only one term in the sum in ( p3| ) contributes to the probability Pmm', 
so that, 



(58) 



V \ / V 



£m-£m'J \\k-k'\5uj 



where we put Cq w £„, e , ~ £m', Vqq' ~ V, \k — k'\ is the "distance" from the 
non-resonant fc'th spin to the resonant fcth spin. 

The total probability of generation of all unwanted states by one pulse in the 
result of non-resonant transitions can be obtained by summation of Pmm' over m' . 
Since each state \m') differs from the state \m) by the flip of one fc' spin, one can 
replace the summation over the states m' to the summation over the spins fc' ^ fc. 
For example, the total probability, fiN-i (here the subscript of fi stands for the 
number of the resonant spin, fc = iV — 1), of generation of unwanted states by the 
initial 7r/2-pulse is, 

(^^) ^^-^ = ^E |^-i-fc^|2 ' ^-[^) ■ 

After the initial 7r/2-pulse (which creates a superposition of two states with equal 
probabilities (Es) from the ground state) the probability of the correct result is 
7'7r/2 = 1 ^ ^J■N-l and the probability of error is PTr/2 = 1 — Ptt/2 = Mjv-i- After 
the first 7r-pulse the probability of error becomes, 

(60) Pi = 1 - 1(1 - f,L-l)il - ^iL-2 -e)- 1(1 - ^iL-l){l - /iL-2). 
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Here we suppose that the values of e = e„ are the same for aU pulses, and Qn — ^ 
for n 7^ 3 and fis = 2ri. 

The probability of error in implementation of the whole logic gate by 7r/2-pulse 
and all 2N — 3 7r-pulses is, 



JV-3 



P = P2N-3 = l--(l-Ai7V-l)(l-MAr-2-e)(l-4AiAr-2-e)(l-M0-£) J^ {^-^^i^^ 



i=l 



N-3 



(61) i(l~Aijv-2)(l-4Mjv-2) n(l-^«) 



where fXi can be derived from (jM) changing A^ — 1 to i. Two last terms in ( p)\ ) 
are connected with two terms in ([45|). The probability of non-resonant transitions 
generated by the third pulse is approximately four times larger (terms with the fac- 
tor 4 in (pi])) than other probabilities because the Rabi frequency for this pulse is 



larger, JI3 = 20 (to keep s the same). Eq. (61 ) takes into account all near-resonant 
transitions, characterized by the parameter e, and non-resonant transitions, char- 
acterized by the parameters fii. There is no need to consider the phases of the 
unwanted states, it is enough to consider only their probabilities. 

Eq. (EW is convenient for choosing optimal parameters for application to the 
scalable solid state quantum computer with a large number of qubits. On the one 
hand, the number of qubits is a scalable parameter in (|6l|), which can be easily 
changed and increased. (For example one can put N — 1000.) On the other 
hand, our approach takes into consideration all significant near-resonant and non- 
resonant processes generating errors in the quantum logic gates. The parameter e 
can be minimized (down to the value e = 0) by the 27rfc-method for any number 
of qubits. However, another parameter, /i, can be significantly decreases only by 
increasing the gradient of the magnetic field, or Slu. Since the number of non- 
resonant transitions is large, their contribution to the probability of errors quickly 
increases with N increasing. As a consequence, for given Su there is a restriction 
on the number of qubits in our computer, N < Nmax, if we want to keep the errors 
below some definite threshold (see section g) . 

8. Computer simulations for finite e 

Here we report our results on computer simulations of a quantum CN gate for 
remote qubits in quantum computer with 200 and 1000 qubits. We have developed 
a numerical code which allows us to study the dynamics of all quantum states with 
the probabilities no less than 10~^ for a spin chain with up to 1000 qubits. All 
frequencies in this and next sections are dimensionless and measured in units of 
J. The values of fl for the results presented in this section are the same for all 
pulses, iln = fl, so that £„ = £(17/2 J) = e for n 7^ 3 and £3 = £(17/4 J). All 
probabilities in this section are doubled. The sum of the probabilities of all these 
states was 2 — O(10~^) (the normalization condition). In this section we suppose 
that £ is large, £ ^ /x (we suppose /x = and Slu = c»), and one can neglect 
the non-resonant transitions, since their contribution to the probability of errors 
is small. The numerical results for the case of finite /i are presented in the next 
section. 

In Fig. 2a, we show the probability of the excited unwanted states after 
implementation of the CiVigg g gate, for N — 200 and fi = 0.14. Because we chose 
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Figure 2. Probabilities of unwanted states. The total number 
of qubits: N = 200; fl = 0.14. (The value of fi is measured in 
units of J.) The number of unwanted states with probabilities 
|C'„p > 10^^ is 7385. The unwanted states are presented in the 
order of their generation, (a) Two "bands" : the upper band and 
the lower band; (b) The upper band of Fig. 2a, magnified; (c) The 
lower band of Fig. 2a, magnified. 



2J, it follows from 



that the closest value of the Rabi frequency, Q,, 
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Figure 3. (a) The ground state of the spin chain; (b-d) Typical 
unwanted states with probabihties P' > 10^''. The horizontal axis 
shows the position of a qubit in the spin chain of A^ = 200 spins. 
The vertical axis shows the ground state, |0), or the excited state, 
|1), of the qubit. Examples of "high energy" (b), "low energy" 
(c) , and "intermediate energy" (d) unwanted states of the whole 
chain. 
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Figure 4. (a) Probabihty, |CoP, as a fimction of Q. (The value 
of Q is measured in units of J.) The total number of qubits, N — 
200; (b) The total number of unwanted states with probabilities 
|C„|2 > 10-*^. 



which satisfies the 27rfc condition is: 




0.1432. 



This value differs slightly from the one used in the simulations whose results are 
presented in Fig. 2a {^ = 0.14) by 3.2 x 10"^. At ft = fit^i there are no the 
near-resonant transitions and the probability of errors is of order of /i, i.e. very 
small. One can see in Fig. 2a that even this small deviation from the 27rfc condition 
results in generating over 7000 unwanted states during the total protocol. On the 
horizontal axis in Fig. 2a, the unwanted states are shown in the order of their 
generation. A total of 7385 unwanted states were generated whose probability 
P' > 10"*^. (In aU figures 2-4 only the states with P' > 10"^ are taken into 
account.) The probability distribution of unwanted states clearly contains two 
"bands". One group of these states has the probability, P' ^ 10^^ (the bold 
"line" near the horizontal axis). The second group of states has the probability 
P' r^ 10~^ (the upper "curve" in Fig. 2a). Fig. 2b shows an enlargement of the 
upper "band" of the Fig. 2a. One can see some sub-structure in this "band" . Fig. 
2c shows the sub-structure in the lower "band" of Fig. 2a. Our simulations show 
the existence of sub-structure in the upper and lower "bands" shown in Fig. 2c. 
There exists a hierarchy in the structure of the distribution function of unwanted 
generated states. Figs. 3(b-d) show the typical structure of unwanted states of 
the spin chain. In Figs. 3(a-d) "0" corresponds to the ground state, |0), of a 
qubit, and "1" corresponds to the excited state, |1), of a qubit. Fig. 3a shows the 
ground state of the spin chain. (All qubits are in their ground state, |0).) P' in 
Fig. 3, is the probability of a state of the whole chain. All unwanted states in Figs 
3(b-d) belong to the upper and lower "bands" shown in Fig. 2c, and they have 
probabilities, P' ~ 10"^. 

Note, that many unwanted states are the high energy states of the spin chain 
(many-spin excitations). Typical unwanted states contain highly correlated spin 
excitations. Fig. 4 shows the probability of the ground state, |CoP, and the total 
number of unwanted states (with probability P' > 10^^) as a function of the Rabi 
frequency, fl. (The point fi = 0.14 corresponds to the results shown in Fig. 2.) 
The maximum value of |CoP and the minimal total number of unwanted states 
correspond to values of fi which satisfy the 2TTk condition for 396 of the total 
number of 7r-pulses, 397. (The third 7r-pulse does not satisfy the 27rfc condition.) 
One can see from Fig. 4 that application of the 27rfc condition significantly improves 
the performance of the quantum Control-Not gate. Fig. 5 shows the ground state 
and the examples of unwanted states for N ~ 1000. In this case, for all 7r-pulses 
n = 0.1, except for the 7r-pulses from 10-th to 40-th for which il = 0.1 + ri, where 
?7 is a random number in the range: —0.05 < i] < 0.05. The value Q = 0.1 



corresponds to the 27rfc condition (44) with k = 10. One can see from Fig. 5 
that the states of the quantum computer are strongly correlated and some of 
them are highly excited. To avoid these unwanted strongly correlated and highly 
excited states it will be necessary to apply both a 27rfc-method and adequate error 
correction codes. 

9. Numerical results for total error 

In the previous section we presented the numerical results when the conditions 
of the 27rfc-method are not satisfied and e ^ /i. In this case one can neglect 
the non-resonant transitions and calculate the error caused by the near-resonant 
transitions by Eq. (|5l|) or (52). One can minimize these errors by choosing £„ = 



for all n (2nk condition). We show in this section that even under the conditions of 
the 27rfc-method the error can be relatively large due to non-resonant transitions, 
and the probability of errors increases when the number of qubits, N, increases. 
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Figure 5. (a) The ground state of the chain; (b-j) Examples of 
unwanted states {N = 1000); O = 0.1 for aU 7r-pulses except for 
the TT-pulses from the 10-th to the 40-th for which il = 0.1 -I- ?7, 
where 77 is a random parameter in the range: —0.05 < 77 < 0.05. 



We present here the results of computer simulations taking into consideration 
both, near-resonant and non-resonant transitions. We test our approximate for- 
mula (pW) by exact numerical solution using Eqs. ( pl| ) and (£2) when the number 
of qubits is not very large {N = 10), so that the total number of states in the 
Hilbert space is 2^" = 1024. As before, the frequencies are dimensionless and 
measured in units of J. The values of e for the results presented in this section 
are the same for all pulses, e„ = £, so that ri„ = il for n ^ 3 and fl^ — 2f7. The 
norm of the wave function is equal (as usual) to unity and 6u> (and fj,) is finite. 

In Figs. 6(a,b) we compare the total probability, P, of generation of unwanted 
states, found from the analytical estimate (pl[), with the result of exact numerical 
solution for the chain containing iV = 10 spins. As follows from these figures, 
there is a very good correspondence between our approximate and exact solutions. 



From Fig. 6a one can see that at e = the value of P has tendency to decrease 
as const / 5lo^ . In Fig. 6b we plot the probability, P, as a function of J7. When 
e ^ ^1, and J <C Sco, the probability of generation of unwanted states is mostly 
defined by the value of e and is almost independent of Scu since e ^ /x. The value 
of Slo fixes the values of the minima in Fig. 6b: for larger Slo the minima in the 
plot in Fig. 6b become deeper. The values of different minima in Fig. 6b indicate 
the contribution of non-resonant processes to the probability P. Since the value 
of Q in Fig. 6b does not change significantly, the contribution of non-resonant 
processes to the probability of errors is approximately the same for all fl. One can 
see that this contribution is negligible in comparison with the contribution of the 
near-resonant processes (defined by e) for all O, except for the small regions of fl, 
where s = 0. 

We should note that one more condition (except for e, /x ^ 1) must be satisfied 
for Eq. (^l|) to be valid. The value of the interaction constant, J, should be small 
in comparison with the difference between the spin frequencies, J <IC Slo. In Figs. 
7a and 7b we plot the probability, P, as a function of Suj for e = and e y^ 0. 
(Since the frequencies are measured in units of J the value 5a; = 1 in Figs. 7(a,b) 
corresponds to Suj ~ J.) One can see that our results are valid only when J <^ Su) 
{5ijj ^ 1 in Figs. 7(a,b)], in spite of the fact that the parameter J /Slo does 
not appear explicitly in (^l|). From Fig. 7b one can see that the probability of 
unwanted states, P, for e ::^ ^ (for large Suo) becomes relatively independent of Slo. 
In this case the value of P is defined by the parameter e which does not depend 
on Slo. 

Equation ( |6l| ) is convenient for choosing the real experimental parameters 
required for operation of the solid state quantum computer. Suppose that we 
are able to correct the errors with the probability less than Pg = 10~^. Our 
perturbation theory allows us to calculate the region of parameters for which the 
probability of error, P, will be less than Pq. In Figs. 8(a-d) we plot the diagrams 
obtained using our perturbative approach for A^ = 10. Inside the hatched areas the 
probability of generation of unwanted states is less than Pp. We should note that 
building the plots like those in Figs. 8(a-d) using exact solution of the problem 
requires a significant number of computer time even for relatively small number 
of qubits {N — 10). On the other hand, the number TV in our perturbation theory 
is the parameter which can be increased without any problems. 

From Figs. 8(c,d) one can see that the sizes of the hatched areas are small, 
much smaller than the areas between the neighboring hatched regions. The total 
size of the hatched regions is mostly defined by the probability, Pq, of tolerant error 
which can be corrected using additional error correction codes. From comparison 
of Fig. 8c for Pq = 10~^ and Fig. 8d for Pq = 10~^ one can see that the sizes of 
the hatched regions increase with Pq increasing. The coordinate of the point A in 
f2 is J7 = rj['^], where the values of k are indicated in the figures. 

In almost all quantum algorithms the phase of the wave function is important. 
We numerically compared the phase of the wave function on the boundaries of the 
hatched regions in Figs. 8(a,b) with the phase in the centers of these regions, 
(where Vt satisfies 27rfc condition, see Eq. (pq), fc = 5 and k = 11) at fixed Soo. 
The deviation in phase is only ^ 0.15%. This is much less that the corresponding 
change in the probabilities of errors (by several orders) . 

In Figs. 9(a,b) we plot the same as in Figs. 8(a,b) but for A^ = 1000. One 
can see that increasing the number of qubits increases the errors and the sizes of 
the hatched areas decrease considerably. From Figs. 8(a-d) and 9(a-b) one can 
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Figure 6. The probability, P, of generation of unwanted states. 
Filled circles are the results of exact numerical solution, dashed 
lines are the analytic estimates (pi]), (a) P as a function of Suj, 
e = 0, f} = Ol^l from {^. (b) P as a function of n, 6oj = 100. 
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Figure 7. The total probability, P, of unwanted states as a func- 
tion of Slo when the value of Slu is comparable with the value of 
the interaction constant, J. (a) fi — fjl^l (e — 0), (b) fi = 0.15 
(e — 0.0039). Filled circles are the results of exact numerical 
solution, solid line is the analytical estimate (|6l|). A^ = 10. 



see that even when the values of fl satisfy 27rfc-method (position of the point A 
in il in Figs. 8(a,b) and extreme left points of the hatched regions in Figs. 8(c,d) 
and 9(a,b)) the error can be relatively large because of the non-resonant processes. 
From comparison of Figs. 8(a,b) for A^ = 10 with Figs. 9(a,b) for N = 1000 one 
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Figure 8. The probability of generation of unwanted states, P, 
at different values of Suj and fi. In the hatched regions P < Pq. 
The different pictures (a) - (d) indicate the different regions in fi 
in different scales, (a), (b), (c) Pq = lO^^ (d) Pq = lO^^. N = 
10. 



can see that the minimal value of Slo required to make the errors small increases 
considerably with TV increasing. 

We now analyze the probability of errors as a function of 6lu. When the value 
of Slo is large enough, the probability of error (and the widths of the hatched 
areas in O) becomes practically independent of Sco. This is because at (5w S> 1 
and at e » ^ the error is mostly defined by s, which is independent of Suj. As 
a consequence, one can, for example, estimate the widths of the hatched areas at 
Suj ^ 1 taking into account only the near- resonant transitions (formula (p2h). 

In order to test our approximate result we solved the problem exactly using the 
eigenstates of the matrices H'^^\ n — 1 , . . . Af , for the parameters which correspond 
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Figure 9. (a) The same as in Fig. 8a but for N = 1000. (b) The 
same as in Fig. 8b but for A^ ~ 1000. The scale in dw shifted to 
the larger values in comparison with Figs. 8a, 8b. 
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Figure 10. The exact solution for the probability of generation 
of unwanted states, P, for the parameters, which correspond to: 
(a) curve AB in Fig. 8a; (b) curve AB in Fig. 8b. The dashed 
lines indicate the solution obtained by the perturbation theory. 
iV = 10. 



to the lower boundary of the hatched regions in the Figs. 8a and 8b (curves AB). 
In Figs. 10a and 10b we compare the total probability of errors, P, obtained using 
the exact and approximate solutions. One can see again, as in Figs. 6a and 6b, 
that there is a good correspondence between the exact and approximate solutions. 
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A similar correspondence can be demonstrated for other points of the parametric 
space {Slli, ft) when the conditions fi ^ J <C (5w are satisfied. 

10. A classical Hamiltonian for quantum computation 

In this section we demonstrate that the process of quantum computation, 
including creation of the entanglement, the dynamics of quantum controlled op- 
erations, and the dynamics of complicated quantum algorithms can be modeled 
using classical Hamiltonians. That is not surprising because the basic quantum 
equations ( |l2[ ) are c-number equations and can be formally considered as an ef- 
fective classical system of equations. The above results show that in some cases, 
effective classical Hamiltonians can be used for simulation of quantum logic opera- 
tions even for large number of qubits. We present here the corresponding effective 
classical Hamiltonian and the Hamiltonian equations of motion in explicit form. 

Formally, we can represent the solution of the Schrodinger equation for the 

Hamiltonian (H)-(Q) in the form: ^(i) = X]ri=o c„(i)|ri), where \n) and En are 
the eigenfunctions and the eigenvalues of the Ising part of the Hamiltonian in (y) 
in the ^-representation, 

(62) Ho\n)^En\n). 
The complex coefficients, c„(i), satisfy the equations, 

2"-l 

(63) icn{t) = ^„c„(t) + Y^ K,fc(t)Cfc(i), 

k=0 

where, 

(64) K,fc = {n\V{t)\k), 

and V{t) is the interaction Hamiltonian in (||). 

Equation (pSh can be written in Hamiltonian form: 

(65) ic„ = — -^, zc„ = -- — , 

dc* dCn 

where Hd is the Hamiltonian of the equivalent "classical" system, 

2"-l 2"-l 

(66) Hcl = ^ En\Cn\'^ + ^ C^Vn,pCp. 

n— n.p—0 

Instead of using complex "classical" variables, c„ and c* , one can introduce the 
following real "classical" variables: a "coordinate", a;„, and a "momentum", p„, 
using, for example, the canonical transformation, 

1 i 

(67) a;„ == — =(c„ + c;), p„ = -— =(c„ - <). 

Using Xn and p„, Eqs. ( pq ) take the familiar classical Hamiltonian form, 

(68) Xn = ^— , Pn = -j^, 
where the equivalent classical Hamiltonian is, 

2"-l -. 2"-l 

(69) Hd{Xn,Pn,t)^ ^ -Y^^n+Pn)+ 2 ^ (^^ ~ iPn)Vn,k(xk + iPk) ■ 

n—0 7i,k—0 
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The Hamiltonian mu) can be written in the form, 



(70) Hcl{Xn,Pn,t) ^ ^ -lJ-{Xn+P' 



2 —1 

^" ' 2 I 2 ^ 



2 

n=0 



n,fc=0 

where we used the relations: Vn^k = VCm ^^^ Re{Vn^k) — Re{Vk^n), Im{Vn,k) = 
^Im(Vk,n). The corresponding classical Hamiltonian equations follow from (pq) 
and ([70|), 



(71) in = EnXn + Y^ [Re{VnM)Pk + 2/m(K,fe) 



Xk\ 



k=0 



2'' 

Pn = -EnPn " ^ [Re{Vn,k)xk " 2Im{Vn^k)Pk]- 
k=0 

For example, the Hamiltonian (|70| ) can be used to calculate the dynamics of the 
quantum Hamiltonian (2). In this case, the matrix elements, Vn,k =/= only for 
those states \n) which differ from the state \k) by a single-spin transitions, and are 
zero otherwise. 

Solutions of Eqs. (|70| ) satisfy the normalization condition, 

2"-l 

(72) Y.i^n{t)+Pl{t)]^l. 

n=0 

Eqs. ( ^0| ) describe the Hamiltonian dynamics of 2^ classical "generalized" one- 
dimensional oscillators. Each of these oscillators is described by two canonically 
conjugate variables, x„ and p„. One can see that 2^ classical harmonic oscillators 
can simulate the behavior of N interacting quantum qubits, including the dy- 
namics of quantum entanglement, complicated quantum logic gates, and quantum 
algorithms. 

Conclusion 

In this paper, we developed an approach, based on small parameters, for the 
simulation of simple quantum logic operations in a nuclear spin quantum computer 
with large number of qubits. We considered a quantum computer which is a one- 
dimensional nuclear spin chain placed in a slightly non-uniform magnetic field, 
oriented in a direction chosen to suppress the dipole interaction between spins. We 
took into consideration the Ising interaction between neighboring qubits. Quantum 
logic operations are implemented by applying resonant electromagnetic pulses to 
the nuclear spin chain. The electromagnetic pulse which is resonant to a selected 
qubit is non-resonant for all other qubits. This raises the problem of minimizing the 
influence of unwanted non-resonant effects in the process of performing a quantum 
protocol. 

We simulate the creation of the long-distance entanglement between remote 
qubits, {N — l)-st and 0-th, in a nuclear spin quantum computer having a large 
number of qubits (up to 1000). We used two essential assumptions: 
1. The nuclear spin chain is prepared initially in its ground state. 
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2. The frequency difference between tlie neighboring spins due to the inhomogene- 
ity of the external magnetic field, Slo, is much larger than the Ising interaction 
constant, J, and J is much larger than the Rabi frequency, fi, i. e. fi <C J "C 5u). 
Using these assumptions, we developed a numerical method which allowed us to 
simulate the dynamics of quantum logic operations taking into consideration all 
quantum states with the probability no less that P = 10^^. For the case when the 
27rfc condition is satisfied (the 7r-pulse for the resonant transition is at the same 
time a 27rfc-pulse for non-resonant transitions), we obtained an analytic solution 
for the evolution of the nuclear spin chain. In the case of small deviations from the 
27rfc condition, the error accumulates and the numerical simulations are required. 
These results are presented in sections ^ and ^. 

The main results of our simulations are the following: 

1. The unwanted states exhibit a band structure in their probability distributions. 
There are two main "bands" in the probability distribution of unwanted states. 
The unwanted states in these "bands" have significantly different probabilities, 
Piow/Pupper "^ 10"'^. Each of these two bands have their own structures. 

2. A typical unwanted state is a state of highly correlated spin excitations. An 
important fact is that the unwanted states with relatively high probability include 
high energy states of the spin chain (many-spin excitations) . 

3. The method developed allowed us to study the generation of unwanted states 
and the probability of the desired states as a function of the distortion of r/ pulses. 
This can be used to formulate the requirements for acceptable errors in quantum 
computation. 

The results of this paper can be used to design experimental implementations 
of quantum logic operations and to estimate (benchmark) the reliability of experi- 
mental quantum computer devices. Our approach can be extended to simulations 
of simple quantum arithmetic operations and fragments of quantum algorithms. 
These simulations are now in progress. 
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